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Abstract. Using realistic hydrodynamical simulations of the 
solar surface convection as 3D, time-dependent, inhomoge- 
neous model atmospheres, the solar photospheric Si abundance 
has been determined to be logesi = 7.51 ± 0.04. This consti- 
tutes a difference of 0.04 dex compared with previous estimates 
based on the ID Holweger-Miiller (1974) model, of which half 
is attributable to the adopted model atmosphere and the remain- 
ing part to the improved quantum mechanical broadening treat- 
ment. As a consequence, all meteoritic abundances should be 
adjusted downwards by the same amount. In particular the me- 
teoritic Fe abundance will be log ep c = 7.46 ± 0.01, in good 
agreement with the recently determined photospheric Fe abun- 
dance (Asplund et al. 2000b). The existing uncertainties un- 
fortunately prevent an observational confirmation of the postu- 
lated effects of elemental migration of metals in the Sun. 
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consistent, 3D, hydrodynamical simulations of the solar sur- 
face convection (e.g. Stein & Nordlund 1989, 1998; Asplund 
et al. 2000a,b, hereafter Paper I and II) and the significance 
of the solar photospheric Si abundance, a re-analysis seemed 
warranted, which is presented here. These convection simula- 
tions have now reached a very high degree of realism. With- 
out relying on any free parameters, such as the mixing length 
parameters, the simulations successfully reproduce the granu- 
lation topology and statistics (Stein & Nordlund 1998), con- 
straints from helioseismology (Rosenthal et al. 1999) and de- 
tailed spectral line shapes, shifts and asymmetries (Paper I and 
II). In particular, line profiles are accurately described without 
invoking any micro- and macroturbulence which are necessary 
in classical ID spectral analyses. As a consequence several pos- 
sible sources of uncertainties (e.g. model atmosphere limita- 
tions, equivalent widths and microturbulence) can be removed 
when deriving elemental abundances, which should result in 
more secure determinations. 



1. Introduction 

To compare the meteoritic and solar photospheric chemical 
compositions it is necessary to have a common reference el- 
ement. The honour is normally given to Si, which is an abun- 
dant element and a natural choice to measure other elemental 
abundances against in meteorites considering the volatility of 
hydrogen. In order to place all other meteoritic abundances on 
an absolute scale to confront with the photospheric abundances, 
one must therefore accurately know the photospheric Si/H ra- 
tio. As a consequence all meteoritic absolute abundances will 
depend on the measured solar photospheric Si abundance. In 
practice, one often utilizes additional elements besides Si to 
anchor the two scales more firmly to each other (e.g. Anders & 
Grevesse 1989). 

Several determinations of the solar photospheric Si abun- 
dance exist in the literature (e.g. Holweger 1973; Lambert & 
Luck 1978; Becker et al. 1980) using the Holweger-Miiller 
semi-empirical model atmosphere (Holweger & Miiller 1974) 
and the use of equivalent widths and the microturbulence con- 
cept. Given the recent progress in the construction of self- 



2. 3D spectral line formation and atomic data 

The procedure is identical to that adopted for a recent deter- 
mination of the solar photospheric Fe abundance (Paper II). 
In particular the model atmosphere is a realistic 3D hydrody- 
namical simulation of the solar surface convection consisting 
of 50 min solar time. For further details of the simulation and 
the 3D spectral line formation the interested reader is referred 
to Paper I and II. 

The Si lines used for the study are the same as in Holweger 
(1973) and are presented in Table [l] together with the necessary 
atomic data. The sample consists of 19 Si I and 2 Sill lines. 
The necessary transition probabilities are taken from Becker et 
al. (1980), which are the same as measured by Garz (1973) but 
corrected using more recent lifetime measurements. We note 
that the lifetime measurements of Smith et al. (1987) in fact 
support the scale of Garz (1973) but until the conflicting life- 
time situation has been clarified, I have chosen to adopt the 
values of Becker et al. (1980) to emphasize the effect of us- 
ing 3D model atmospheres; new improved measurements of 
the branching ratios and lifetimes of the Si lines seem jus- 
tified. The collisional broadening data for both Si I and Si II 
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Table 1. The adopted line data for the Si I and Si II lines 



Species 


Wavelength 


XI 






log esi 




[nm] 


[eV] 




[pm] 




Si I 


564.56130 


4.93 


-2.04 


3.40 


7.52 


Si I 


566.55550 


4.92 


-1.94 


4.00 


7.45 


Si I 


568.44840 


4.95 


-1.55 


6.00 


7.47 


Si I 


569.04250 


4.93 


-1.77 


5.20 


7.51 


Si I 


570.11040 


4.93 


-1.95 


3.80 


7.48 


Si I 


570.84000 


4.95 


-1.37 


7.80 


7.48 


Si I 


577.21460 


5.08 


-1.65 


5.40 


7.51 


Si I 


578.03840 


4.92 


-2.25 


2.60 


7.52 b 


Si I 


579.30730 


4.93 


-1.96 


4.40 


7.55 


Si I 


579.78560 


4.95 


-1.95 


4.00 


7.55 


Si I 


594.85410 


5.08 


-1.13 


8.60 


7.48 


Si I 


674.16280 


5.98 


-1.65 


1.60 


7.55 


Si I 


697.65130 


5.95 


-1.07 


4.30 


7.56 b 


Si I 


703.49010 


5.87 


-0.78 


6.70 


7.54 b 


Si I 


722.62080 


5.61 


-1.41 


3.60 


7.51 


Si I 


768.02660 


5.86 


-0.59 


9.80 


7.54 


Si I 


791.83840 


5.95 


-0.51 


9.50 


7.50 b 


Si I 


793.23480 


5.96 


-0.37 


9.70 


7.53 b 


Si I 


797.03070 


5.96 


-1.37 


3.20 


7.59 


Sill 


634.71090 


8.12 


0.30 


5.60 


7.48 


Sill 


637.13710 


8.12 


-0.00 


3.60 


7.43 



1.2 
1 .0 



From Holweger (1973). Note that W\ is only listed here to allow 
easy identification in Fig. ^ and is not used for the abundance 
determinations 

Lines which are entered with half weight into the final abundance 
estimate due to uncertainties in broadening, oscillator strengths or 
suspected blends 



lines have been computed following the recipes of Anstee & 
O'Mara (1991, 1995), Barklem & O'Mara (1997), Barklem et 
al. (1998) and Barklem (2000, private communication). If one 
instead would embrace the classical collisional treatment of 
Unsold (1955) without any damping enhancement factors the 
final Si abundance would be 0.02 dex higher with a slightly in- 
creased scatter. Stark broadening and radiative broadening was 
included with data taken from the VALD database (Kupka et al. 
1999), although Stark broadening is only of some significance 
for Si I 730.49 and 768.03 nm. The Si abundance has been de- 
termined from profile fitting of the individual Si lines through 
a x 2 -analysis in a similar fashion to the study of Nissen et al. 
(2000) rather than equivalent widths. The solar atlas of Brault 
& Neckel (1987) and Neckel (1999) provided the observed 
spectrum. Fig. [TJ shows an example of the achieved agreement 
between predicted and observed line profiles. It should be em- 
phasized that this has been accomplished without invoking any 
free parameters besides the elemental abundance; the result 
is clearly superior to what is possible using ID model atmo- 
spheres such as the Holweger-Muller (1974) model. 



0.2 
0.0 




A = 768.027 [nm] 
X = 5.863 [eV] 
log gf = -0.590 
log e = 7.540 



768.00 768.02 768.04 

Wavelength [nm] 



Fig. 1. The spatially and temporally averaged (diamonds) Si I 
768.0 nm line calculated using the 3D inhomogeneous model 
atmosphere compared with the observed profile (solid). The 
predicted profile has been convolved with a sinc-function to ac- 
count for the finite spectral resolution of the solar atlas (Paper 
I and II). Since no free parameters but the Si abundance enter 
the spectral synthesis, the good agreement supports the con- 
clusion that the macroturbulence concept required in classical 
ID analyses is not necessary when properly accounting for the 
convective Doppler shifts in the hydrodynamical simulations 

3. The solar photospheric Si abundance 

Fig. |2] reveals no clear trend in the derived individual Si abun- 
dances with line strength; the span in excitation potential is un- 
fortunately too small to clearly delineate any possible trend. 
When including all Si I and Si II lines with equal weight, the 
estimated solar photospheric Si abundance is log egi = 7.51 ± 
0.040, where the quoted uncertainty is the standard deviation. 
If equivalent widths from the literature (Holweger 1973) had 
been used instead of profile fitting, the mean abundance would 
have remained the same but the scatter would have increased 
somewhat. When comparing line profiles it was noticed that a 
few lines are less suitable for abundance determinations. The 
most likely reasons for the poor agreement in the profile fitting 
appear to be suspected blends or erroneous broadening treat- 
ment. Such lines enter the final abundance estimate with half 
weight as marked in Table [l]. 

The re-analysis of the solar Si lines then results in 

logesi = 7.51 ±0.04, 

which is only slightly less than the value (7.55 ± 0.05) pub- 
lished by Becker et al. (1980) with the same set of transi- 
tion probabilities. Half of this difference is due to the adop- 
tion of a 3D hydrodynamical model atmosphere instead of the 
Holweger-Muller (1974) model atmosphere and the remaining 
part to the improved quantum mechanical broadening treatment 
instead of Unsold (1955) pressure broadening. Clearly the Si 
lines are little sensitive to the detailed atmospheric structure 

1 On the customary logarithmic abundance scale defined to have 
loge H = 12.00 
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Fig. 2. The individual Si abundances derived from the profiles 
of Si I (solid circles) and Si II (open circles) lines as a function 
of the equivalent widths (taken from Holweger 1973); it should 
be emphasized however that equivalent widths have not been 
used for the abundance determinations. The dashed line rep- 
resents a least square fit to the Si abundances. The lack of an 
apparent trend supports the conclusion that no microturbulent 
velocity is necessary to include in the spectral synthesis when 
correctly including the Doppler shifts inherent in the convec- 
tion simulations 



and emphasizes that the main uncertainty is likely to be the 
absolute gf '-scale, although the small scatter suggests that the 
relative measurements by Garz (1973) are of good internal con- 
sistency. Departures from LTE will probably not significantly 
influence the derived mean abundance, considering the rela- 
tively small scatter for the Si I lines and since ionization bal- 
ance seems to be essentially fulfilled; some caution in this mat- 
ter should still be exercised however, since only two Si II lines 
enter the analysis. Additional Si II lines would certainly be very 
helpful in securing the results further. 

4. The photospheric and meteoritic abundance scales 

A re-analysis of the solar photospheric Si abundance using 3D, 
time-dependent, inhomogeneous model atmospheres have re- 
vealed only a difference of 0.04dex (10%) with previously 
published estimates. As described in Sect, [jjthis also suggests 
that the meteoritic abundances remain largely unaltered com- 
pared with the commonly adopted values (Grevesse & Sauval 
1998). If only Si is used to place the meteoritic and photo- 
spheric abundances on a common scale, the meteoritic Fe abun- 
dance for example will be log ep e = 7.46 ± 0.01 (an additional 
0.01 dex decrease may come from the existing difference in 
adopted meteoritic and photospheric Si abundances, Grevesse 
& Sauval 1998). It should be noted though that several other el- 
ements are often utilized as well (Na, Mg, Si, Ca, V, Cr, Co, Ni, 
Y, Zr, Nb and Mo, Anders & Grevesse 1989) for the purpose 
and it is not clear whether all these elements will behave sim- 
ilarly to Si when changing from the Holweger-Miiller (1974) 



model atmosphere to the 3D solar convection simulations form- 
ing the basis of this work. However, given the similarity in this 
respect between Si and Fe (Paper II) it seems reasonable that all 
meteoritic abundances indeed need to be adjusted downwards 
by about 0.04 dex. 

The difference between the meteoritic Fe abundance ar- 
rived at here and the recently determined photospheric abun- 
dance (Paper II) is only 0.02 dex for Fe I (log e Fc = 7.44±0.05) 
and 0.01 dex for Fell (logeFo = 7.45 ± 0.10) lines, when 
solely relying on the Si results presented here for the absolute 
scale for the meteoritic abundances. It is worthwhile to keep in 
mind, however, that the photospheric Fe results using 3D model 
atmospheres are based on the assumption of LTE. For some 
ID model atmospheres there are indications that low-excitation 
Fe I lines may be influenced by departures from LTE at the level 
of Aloge F e k 0.1 dex (Solanki & Steenbock 1988; Shchukina 
2000, private communication), which would, if assumed to be 
valid also in the 3D case, degrade the agreement between the 
3D LTE results for Fe I and Fe II (Paper II). However, until de- 
tailed 3D NLTE calculations have been performed such a con- 
clusion appears to be premature since the NLTE results are sen- 
sitive to the adopted model atmosphere. Rather the consistency 
between the meteoritic value and the abundances derived from 
both weak and strong Fe I and Fe II lines and the absence of 
any trend in abundances with excitation potential may be inter- 
preted as NLTE effects being insignificant; in any case the Fe II 
results should be robust against departures from LTE (Solanki 
& Steenbock 1988; Shchukina 2000, private communication). 

The difference between the here estimated meteoritic and 
photospheric Fe abundances is slightly smaller than the value 
of 0.04 dex (10%) expected due to diffusion and gravitational 
settling (e.g. Vauclair 1998), similarly to the situation for 
He which has been uncovered by helioseismology (e.g. Basu 
1998). However, the absolute transition probabilities for Fe and 
Si are unfortunately still not of sufficient accuracy to definitely 
reveal such a small effect, even though the atmospheric mod- 
elling now has answered the challenge. Furthermore, and more 
importantly, the absolute abundances in meteorites have been 
obtained by enforcing the average differences between the pho- 
tospheric and meteoritic abundances for some 10 metals to van- 
ish. It is therefore not surprising that the average difference 
in abundances for the 59 elements (excluding e.g. Li, Be, B, 
which may have been depleted) is < 0.01 dex (Grevesse & 
Sauval 1998). Thus, in order to observationally prove elemental 
migration for the Sun through spectroscopy, differential diffu- 
sion effects must be identified, which are expected to be signif- 
icantly less than the overall depletion of metals by about 10%. 
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